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Chapter  1 


INTRODUCTION 


1.1  Introduction 

Ordinary  Portland  Cement  (OPC)  is  a  widely  used  construction  material.  In 
airbase  facilities,  it  is  used  to  construct  the  runways  and  parking  aprons.  During  the 
takeoff  and  landing  of  aircrafts,  the  high  temperature  exhaust  gas  splashed  to  the  concrete 
pavement  from  the  modem  vectored  thmst  engines  (VTE)  or  Auxiliary  Power  Unit 
(APU).  The  APU  is  a  low-power  gas  turbine  that  provides  compressed  air,  from  a  load 
driven  compressor,  for  starting  the  main  engines  and  for  operating  auxiliary  systems 
during  ground  maintenance.  In  particular,  the  exhaust  gas  temperature  of  a  VTE  could 
rapidly  reach  over  700  °C,  and  the  corresponding  exhaust  velocities  could  go  beyond 
1800  ft/s.  The  pavement  is  subjected  to  extremely  rapid  transient  high  temperature 
loadings  as  well  as  thermal  cycles  of  heating  and  cooling.  Figure  1.1  shows  a  typical 
damaged  area  in  an  F/A-18  parking  apron.  The  origin  of  the  problem  is  the  F/A-18 
auxiliary  power  unit,  which  is  located  at  the  bottom  of  the  fuselage  near  the  tail  section. 
The  exhaust  gas  temperature  isotherms  on  the  pavement  top  are  almost  circular  and  show 
a  decrease  along  the  radius. 

Internal  microcracks  nucleate  and  propagate  in  concrete  due  to  these  severe 
thermomechanical  loadings.  The  damage  significantly  affects  the  residual  effective 
moduli  and  strength  of  the  pavement  concrete.  Situations  are  much  worse  to  airfield 
concrete  pavements  for  the  future  vertical  take-off/landing  (VSTOL)  aircraft  (i.e.  the 
Joint  Strike  Fighter)  with  vectored  thrust  engines  (VTE)  [Hooker  (1981)].  The  damage  of 
airfield  concrete  pavements  caused  by  the  aircraft  exhaust  is  an  expensive  problem  that  is 
expected  to  get  worse.  The  damage  occurs  in  the  upper  1/4  to  1/2  inch  of  the  surface  and 
results  in  debris  that  can  severely  damage  the  aircraft  engines.  The  loose  cement  paste 
and  aggregate  may  be  easily  ingested. 


Damage  of  airfield  concrete  pavement  could  be  in  the  form  of  thermomechanical 
spalling,  delamination,  and  splitting.  Spalling  refers  to  small,  thin  pancake-shaped 
pieces  that  explode  into  the  air  because  of  the  very  large  steam  pressure  in  pores,  high 
thermal  gradients  and  large  compressive  thermal  stresses  (parallel  to  the  surface  in  the 
pavement).  Specifically,  compressive  stresses  produce  not  only  concrete  crushing,  but 
also  a  bulging  instability  (similar  to  buckling)  of  the  top  layer.  Under  these  severe 
thermomechanical  loadings,  internal  microcracks  nucleate  and  propagate  in  concrete, 
leading  to  progressive  delamination  at  some  shallow  depth  beneath  the  pavement  surface. 
Figure  1.2  is  the  explosive  spalling  occurring  in  saturated  concrete  heated  to  simulate 
VSTOL  jet  exhaust  (Kodres,  1996).  During  the  heating  procedure  concrete  pavement  is 
under  high  level  compressive  thermal  stress  loading.  The  microcracks  may  propagate  in 
the  plane  parallel  to  the  pavement  surface.  This  mechanism  will  cause  another  type  of 
delamination.  However,  during  the  cooling  procedure,  the  dry-shrinkage  strain  can  not 
recover.  A  high  level  tensile  thermal  stress  will  be  generated.  In  this  case,  the 
microcracks  will  propagate  perpendicular  to  the  pavement  surface.  Finally,  the  pavement 
will  be  splitted.  Since  concrete  is  weak  under  tension,  it  is  important  to  investigate  the 
behavior  of  concrete  under  tensile  loadings.  If  the  pavement  is  delaminated  in  horizontal 
plane  and  splitted  in  vertical  direction  relative  to  the  pavement  surface,  thin  plate-like 
pieces  or  lamina  that  flake  or  peel  off  from  the  damaged  concrete  surfaces.  This 
phenomenon  is  referred  to  as  scaling.  Figure  2.3  is  the  pavement  scaling  due  to  F/A-18 
exhaust  at  the  Naval  Air  Station,  North  Island,  San  Deigo  (Kodres,  1996). 


1.2  Composition  of  concrete 

In  its  most  general  definition,  concrete  is  a  composite  material  composed  of 
coarse  granular  material  (the  aggregate  or  filler)  embedded  in  a  hard  matrix  material  (the 
cement  or  binder)  that  fills  the  space  between  the  aggregates  and  binds  them  together 
(Mindess  and  Young,  1981).  Concrete  is  popular  because  it  is  economical,  and  features 
desirable  casting  and  in-situ  forming  capabilities.  However,  hardened  concrete  is  both 
chemically  and  mechanically  complex.  It  has  been  suggested  that  at  least  three  main 
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structural  levels  or  scales  exist  for  concrete,  and  that  each  of  these  levels  exhibits  unique 
characteristics  with  varying  degrees  of  complexities  (Wittmann,  1983) 

The  structure  and  characteristics  of  the  hardened  cement  paste  are  the  main 
factors  considered  at  the  micro-level.  When  mixed  with  water,  the  constituent 
compounds  of  Portland  cement  undergo  a  series  of  complex  chemical  reactions 
(hydration)  which  lead  to  the  eventual  hardening  of  the  concrete.  These  reactions  are 
always  exothermic  (i.e.,  they  generate  heat).  The  amount  of  heat  generated  during  these 
reactions  results  in  thermal  stresses  which  usually  lead  to  cracking  within  the  matrix. 
After  the  primary  hydration  process  is  complete,  the  new  solids  formed  are  called  the 
hydration  products. 

The  characteristic  features  of  the  heterogeneous  structure  of  the  meso-level 
consist  of  inclusions  (hard  aggregate),  pores,  and  pre-existing  cracks  embedded  in  the 
hardened  cement  paste  (matrix)  (Figure  1.4).  In  normal  weight  concrete,  randomly 
dispersed  aggregates  of  varying  sizes  and  shapes  are  tougher  than  the  cement  matrix  and, 
consequently,  act  as  crack  arrestors.  It  is  the  geometry  of  the  mesostructure  (the 
dispersion  and  volume  fraction  of  constituent  phases,  the  specific  area  of  weak  interfaces, 
and  the  size  and  distribution  of  pre-existing  defects)  which  affects  the  overall  mechanical 
response  and  ultimately  determines  the  rupture  strength  of  the  material.  A  well- 
documented  mesoscale  characteristic  is  that  the  aggregate-cement  paste  interface  is  often 
the  weakest-link  in  the  concrete  composite.  Therefore,  the  cracks  formed  at  this  interface 
have  the  most  significant  influence  on  the  material  behavior.  This  influence  is 
macroscopically  observed  as  the  nonlinearity  in  the  stress-strain  curve. 

At  the  macro-level,  the  size  of  the  specimen  is  considered  to  be  large  enough 
(i.e.,  possessing  a  large  number  of  heterogeneities)  so  that  the  material  can  be  considered 
statistically  homogeneous;  i.e.,  stress/strain  fluctuations  are  not  significant. 
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1.3  Literature  review 


An  improved  understanding  of  the  underlying  deformation  and  damage 
mechanisms  will  help  the  pavement  design  and  prolong  the  usable  life  of  airfield 
pavements.  Research  on  micromechanical  damage  level  is  needed  to  better  understand 
the  basic  physical  phenomena  responsible  for  these  changes. 

High-temperature  thermomechanical  behavior  and  properties  of  concretes  have 
been  studied  in  the  past  thirty-five  years  primarily  due  to  the  interests  in  heat  resistance 
of  concretes  in  fire  engineering,  conventional  and  nuclear  electric  power  plants.  For  an 
extensive  literature  review  on  high-temperature  behavior  of  Portland  cement  and 
refractory  concretes,  we  refer  to  Bazant  et  al.  (1982).  Some  important  aspects  of 
temperature-dependent  properties  of  OPC  concretes  at  elevated  temperatures  discussed 
by  Bazant  et  al.  (1982)  include  dehydration  of  cement,  increase  of  porosity,  moisture 
content,  thermal  expansion  and  shrinkage,  pore  steam  pressure,  loss  of  strength,  thermal 
cracking  due  to  thermal  incompatibility,  degradation  of  elastic  moduli,  thermal  creep, 
heat  capacity,  thermal  conductivity  and  thermal  diffusivity,  and  explosive  thermal 
spalling  due  to  excessive  pore  steam  pressure,  etc. 

The  two  principal  culprits  for  the  damage  of  concrete  pavement  are  the 
microcracks  and  the  high  pore  pressure.  Accordingly,  there  are  two  types  of  thermo¬ 
micromechanical  damage  models  of  concrete.  One  type  of  damage  model  emphasizes  the 
development  of  the  microcracks  in  the  concrete.  Krajcinovic  and  Fanella  (1986)  proposed 
a  constitutive  theory  for  concrete  based  on  the  geometry  of  its  mesostructure  and  the 
actual  kinetics  of  the  microcrack  evolution.  Fanella  and  Krajcinovic  (1988)  develop  a 
general  three-dimensional  micromechanical  constitutive  theory  for  plain  concrete 
subjected  to  uniaxial  and  triaxial  compressive  loads.  Ju  and  Lee  (1991)  present  three- 
dimensional  micro-mechanical  anisotropic  damage '  models  for  microcrack-weakened 
brittle  solids.  Lee  and  Ju  (1991)  developed  three-dimensional  self-consistent  damage 
models  for  brittle  solids  under  compressive  triaxial  loadings  based  on  micromechanics 
and  microcrack  geometry  within  a  representative  volume.  Ju  and  Zhang  (1998)  proposed 
a  simple  one-dimensional  thermo-micromechancal  constitutive  and  damage  model  for 
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airfield  concrete  pavement  under  rapid  transient  high  temperature  loadings  due  to 
vectored  thrust.  Chandra  and  Krauthammer  (1996)  present  a  micromechanical  damage 
model  for  a  brittle  solid  capable  of  taking  into  account  the  effect  of  high  loading  rate  (in 
the  form  of  stress-rate). 

By  experimental  program,  Zoldners  (1960,  1971),  Zoldners  et  al.  (1963),  and 
Zoldners  and  Wilson  (1973)  studied  effects  of  sustained  and  cyclic  high  temperatures  on 
concretes  containing  different  aggregates.  Khoury,  Grainger  and  Sullivan  (1985a,  1985b, 
1986)  conducted  research  on  transient  (free  and  load-induced)  thermal  strain  of  concrete, 
including  first  heating  to  600  °C  and  first  cooling  from  600  °C  under  load.  Dias,  Khoury 
and  Sullivan  (1990)  discussed  mechanical  properties  of  hardened  cement  paste  exposed 
to  temperatures  up  to  700  ”C.  Moreover,  Baluch  et  al.  (1989)  studied  the  degradation 
(microcracking)  of  concrete  bridge  decks  in  Saudi  Arabia  due  to  incompatible 
coefficients  of  thermal  expansion  between  aggregates  and  cement  paste.  Although  the 
aforementioned  studies  (on  high-temperature  behavior  of  concrete)  in  the  literature  are 
quite  useful,  they  concentrate  on  relatively  slow  heating  rates  (e.g.,  less  than  10  ®C/min.) 
of  thermomechanical  properties  of  concretes.  By  contrast,  our  objective  here  is  to  conduct 
preliminary  investigation  of  the  transient  high-temperature  (say  700  °C),  fast  heating  rate 
(say  500  °C/min.)  thermomechanical  behavior  of  ordinary  Portland  cement  concrete 
airfield  pavements  with  various  mixtures  of  cement  pastes  and  aggregates. 

Other  scholars  pay  attention  to  the  high-level  pore  pressure  in  the  concrete 
pavement.  Concrete  is  a  material  that  contains  extremely  fine  pores  with  a  great  amount 
of  water.  Therefore,  heating  of  concrete  will  produce  significant  pore  pressure  that  may 
induce  transfer  of  moisture  through  concrete  and  eventual  drying.  The  high  pore  pressure 
and  the  migration  of  moisture  and  exhaust  caused  by  heating  of  concrete  pavement  due  to 
the  exhaust  gas  from  the  aircraft  is  of  considerable  concern  for  predicting  the  damage  of 
airfield  concrete  pavement.  There  exist  some  studies  on  the  mass  transport  and  heat 
transfer  within  the  stagnation  zone.  For  instance,  Stevens,  Pan  and  Webb  (1992),  and 
Pan,  Stevens  and  Webb  (1992)  experimentally  studied  the  turbulent  flow  structure,  local 
heat  transfer  and  effects  of  nozzle  configuration  on  transport  in  the  stagnation  zone  of 
axisymmetric,  impinging  free-surface  liquid  jets.  In  addition,  Rish  and  McVay  (1994) 
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investigated  the  moisture  release  and  migration  in  concrete  due  to  vectored  thrust. 
Specifically,  Rish  and  McVay  (1994)  conducted  research  on  mass  transport  under 
isothermal  conditions,  and  measured  pore  pressures  during  transient  heating. 

There  are  several  different  approaches  to  this  analysis.  One  approach  is  to  assume 
the  mass  flow  of  the  liquid  phase  is  negligible  and  treat  the  problem  as  a  flow  of  gas 
through  a  porous  concrete.  Sahota  and  Pagni  (1979)  used  this  approach  to  examine  the 
thermodynamics  of  moisture  in  the  pores  of  a  concrete  structure  subjected  to  fire.  They 
solved  the  conservation  equations  for  energy  and  species,  employing  a  finite  difference 
scheme.  Permeability  was  used  as  a  parameter.  However,  they  stopped  the  parametric 
study  when  pore  pressures  reached  10  atm.  Pressures  above  10  atm  “are  not  expected  due 
to  microcracking  and  cracking”.  An  empirical  approach  was  adopted  by  Bazant  and 
Thonguthal  (1978)  to  study  the  heating  of  concrete.  By  fitting  data,  they  developed 
relationships  between  pore  saturation  and  pore  pressure  and  between  permeability  and 
pore  pressure.  These  equations  were  coupled  with  conservation  of  energy  and  mass,  and 
solved  using  finite  element  method.  An  unsealed  concrete  slab  was  modeled  (Bazant  and 
Thonguthal,  1979),  with  surface  temperature  increasing  at  the  rate  of  80  "C/min.  At  this 
rate,  the  surface  temperature  reached  150  °C  in  slightly  less  than  2  minutes,  similar  to  the 
rate  at  which  an  APU  heats  an  airfield  pavement.  The  predicted  maximum  pore  pressure 
was  about  5  atm. 

Another  approach  to  calculating  moisture  migration  through  high-resistance 
porous  media  is  to  use  the  methods  developed  to  study  the  migration  of  moisture  through 
permeable  media  such  as  soil.  The  flow  of  both  liquid  and  gaseous  phases  is  accounted 
for.  Applying  theoretical  correlations  between  capillary  pressure  and  pore  saturation,  and 
relative  permeability  and  pore  saturation  (Van  Genuchten,  1980).  Doughty  and  Pruess 
(1990)  employed  this  method  to  determine  the  feasibility  of  constructing  a  geologic 
repository  for  high-level  nuclear  waste.  Kodres  (1996)  used  a  constitutive  relationship  to 
characterize  the  flow  of  high-velocity,  compressible,  heated  gasses  through  a  high- 
resistance  porous  medium.  The  resistance  to  flow  was  modeled  with  friction  coefficients 
analogous  to  the  method  used  for  simpler  geometries.  Pore  pressure  generated  by  water 
and  air  in  the  pores  of  airfield  concrete  being  heated  by  a  jet  exhaust  was  predicted. 
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Abdel-Rahman  and  Ahmed  (1996)  predicted  the  heat  and  mass  transfer  in  concrete  wall 
exposed  to  fire.  Dayan  and  Gluekler  (1982)  treated  the  concrete  as  a  porous  material 
containing  water,  vapor  and  air.  They  studied  the  heat  and  mass  transfer  within  a  surface 
of  heated  concrete  slab.  The  most  cited  experimental  study  is  still  the  works  of  England 
and  Ross  (1970)  and  England  and  Sharp  (1971).  They  heated  one  end  of  a  0.6  m 
cylindrical  concrete  specimen  and  measured  pressure  and,  using  resistance  meters, 
moisture  along  its  length.  The  sides  of  the  specimen  were  insulated  to  ensure  one 
dimensional  heat  and  mass  transfer.  The  hot  end  reached  150  ®C.  Although  time  and 
space  gradients  of  temperature  were  very  small,  this  work  is  of  special  interest  because  it 
provides,  in  part,  an  experimental  confirmation  of  the  hypothesis  behind  the  airfield 
erosion  analysis. 

I  Chemical  degradation  of  concrete  due  to  hydraulic  fluid,  lubricating  oils  also 

contribute  to  the  damage  of  concrete  pavement.  The  U.S.  Department  of  Defense  has 
seen  an  increase  in  airfield  concrete  apron  distress  in  the  form  of  surface  scaling  when  the 
aprons  are  exposed  to  cyclic  heat,  spilled  lubricants,  and/or  hydraulic  fluids  (McVay, 
Smithson  and  Manzione,  1993).  Chemical  analysis  of  the  damaged  concrete  reveals  that 
the  spilled  fluid  from  the  engines  are  undergoing  hydrolysis  (breakdown)  accompanied 

I  by  the  consumption  of  calcium  hydroxide,  plus  hydrated  silicate  and  aluminate  phases.  A 

I  thermo-chemical-micromechanical  damage  model  is  to  be  derived  as  the  future  research. 


1.4  Scope 

In  Chapter  2,  a  thermo-micromechanical  damage  model  is  proposed  for  the 
airfield  concrete  pavement  under  extremely  rapid  heating  and  cooling  processes  due  to 
high  temperature  exhaust  gas  from  the  vectored  thrust  engines.  This  is  typical  of 
advanced  aircraft  during  their  short  take-off  and  vertical  landing  routines.  Since  the 
temperature  range  in  this  research  is  quite  large,  and  thermal  properties  of  concrete  are 
usually  varying  with  the  change  of  temperature,  thermal  properties  are  treated  as 
functions  of  temperatures.  The  temperature  distribution  of  the  airfield  concrete  pavement 
is  assumed  to  be  independent  of  the  stress  field,  and  is  numerically  calculated  by  the 
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explicit  finite  difference  method,  We  consider  the  purely  elastic  strain,  the  stable  crack- 
induced  strain,  the  unstable  crack-induced  strain,  the  void-induced  strain,  the  vapor 
pressure-induced  strain,  the  thermal  expansion  strain,  and  the  dry-shrinkage  strain.  Based 
on  the  strain  analysis,  a  three-dimensional  constitutive  equation  is  derived  for  concrete 
under  high  temperature  loading.  In  addition  to  the  constitutive  relations,  three  types  of 
failure  mechanisms  are  discussed.  The  pavement  failure  mainly  consists  of  three 
mechanisms:  the  delamination  due  to  horizontal  void-induced  cracks;  the  delamination 
due  to  horizontal  interfacial  cracks;  and  the  splitting  due  to  vertical  cracks.  This  damage 
occurs  in  the  upper  surface  and  results  in  debris  that  can  severely  damage  the  aircraft 
engines.  Within  the  framework  of  linear  thermo-elasticity,  the  stress  distribution  as  a 
function  of  location  and  time  is  also  calculated  by  the  finite  difference  method. 
Furthermore,  based  on  the  ASME  Steam  Tables  and  the  temperature  histories,  we  are 
capable  of  predicting  the  pore  pressure  fields  and  histories.  For  demonstration,  a  general 
one-dimensional  problem  is  numerically  simulated.  Both  normal-weight  aggregate 
concrete  and  light-weight  aggregate  concrete  are  considered.  The  parametric  study  of  the 
concrete  can  lead  to  the  optimum  design  of  the  pavement. 

In  Chapter  3,  the  coupled  governing  equations  of  heat  and  mass  transfer  in 
concrete  pavement  are  derived.  Concrete  is  a  material  that  contains  extremely  fine  pores 
with  water.  Therefore,  heating  of  concrete  will  produce  significant  pore  pressure  that  may 
induce  the  transport  of  moisture  through  concrete  and  the  eventual  drying.  Meanwhile, 
the  migration  of  moisture  through  concrete  may  appreciably  contribute  to  the  heat 
transfer  and  affect  the  temperature  distribution  inside  the  concrete.  The  pore  pressure 
depends  on  the  pore  space  available  to  water,  which  is  known  to  change  in  response  to 
heating  and  other  factors.  The  material  properties  involved  in  the  diffusion  problem  of 
moisture  and  heat  transfer  are  strongly  variable.  Therefore,  analysis  of  pore  pressure, 
moisture  content  and  temperature  in  heated  concrete  necessitates  a  numerical  approach. 
The  objective  of  this  study  is  to  develop  a  finite  difference  program  to  accomplish  this 
task.  For  the  case  of  jet  exhaust  impinging  on  the  concrete  pavement,  migration  of 
exhaust  through  the  concrete  is  also  involved  besides  the  transfer  of  heat  and  moisture. 
The  high  pore  pressure  and  the  migration  of  moisture  and  exhaust  caused  by  heating  of 
concrete  pavement  due  to  the  exhaust  gas  from  the  aircraft  is  of  considerable  concern  for 
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predicting  the  damage  of  airfield  concrete  pavement.  To  model  the  mass  (moisture  and 
exhaust)  and  heat  transfer  through  concrete,  coupled  partial  differential  equations 
concerning  the  conservation  of  moisture,  exhaust  and  energy  must  be  solved.  The  effects 
of  Soret  flux  to  the  moisture  flux  and  gas  flux  are  considered.  Water  vapor  pressure  and 
gas  pressure  are  chosen  to  be  the  driving  forces  that  govern  moisture  transfer  and  gas 
transfer,  in  order  to  eliminate  the  effects  of  temperature  gradient  to  the  moisture  flux  and 
gas  flux.  Dehydration  of  water  in  concrete  at  high  temperature  is  also  considered. 
Empirical  relations  among  water  vapor  pressure,  temperature  and  water  concentration  in 
the  concrete  are  employed.  The  moisture  and  exhaust  gas  transfer  in  concrete  are  also 
investigated.  Semi-implicit  finite  difference  scheme  is  applied  to  calculate  temperature 
and  pore  pressure  distributions.  The  stability  criteria  of  the  numerical  simulation  is  also 
derived. 
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